A Langevin canonical approach to the dynamics of two level systems. I. 

Populations and coherences. 
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A canonical framework for two-level systems coupled to a bath of harmonic oscillators is developed 
to extract, from a stochastic dynamics, the thermodynamic equilibrium values of both the population 
difference and coherences. The incoherent and coherent tunneling regimes are analyzed for an Ohmic 
environment from moderate-to-high temperatures. The corresponding numerical results issued from 
solving a non-linear coupled system are fitted to approximate path-integral analytical expressions 
in order to determine the different time scales governing both regimes. 
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I. INTRODUCTION 

The description of many phenomena in terms of a 
two level system (TLS) can be found in different ar- 
eas of chemical physics (chiral molecules l], electron 
transfer reactions 0, H[), quantum optics 5[, high energy 
physics @ and quantum computation , among many 
others fields of research. However, the isolated TLS is not 
complex enough to provide a richer and more complete 
image of the basic microscopic processes underlying its 
dynamics. For example, very often the system is coupled 
to a more extended system or environment consisting of 
many degrees of freedom usually represented by an infi- 
nite set of harmonic oscillators. Thus, several extensions 
of the isolated TLS have been taken into account along 
the years to account properly for dissipative and stochas- 
tic processes due to interactions with an environment 
at finite temperature. In general, the TLS is linearly 
coupled to the coordinates of a bath of nonintcracting 
oscillators, whose properties are encoded in their spec- 
tral density Q (the spin-boson description). Several ap- 
proaches within the matrix density formalism have been 
proposed and implemented to determine the time evolu- 
tion of the dissipative TLS such as path-integral methods 
or the Bloch Redficld formalism 0-Qji. In the path- 
integral formalism, the so-called nonintcracting blip ap- 
proximation (NIBA) has been shown to be a good ap- 
proximation especially in the regimes of strong damping 
and/or high temperatures providing analytical results to 
the incoherent and coherent tunneling rates Q. Varia- 
tional calculations have also been carried out for both the 
symmetric (or unbiased) [l6| and asymmetric (or biased) 
cases fvj\ . Alternatively, following an approach to the 
dynamics of the TLS due to Feynman in his elegant dy- 
namical theory of the Josephson effect [l8| , it was shown 
that classical and quantum mechanics may be embedded 
in the same Hamiltonian formulation by using complex 
canonical coordinates fl9l . [2pj . In particular, for elec- 
tronically non-adiabatic processes, Meyer and Miller 1211 
(following some earlier work by McCurdy and Miller [22[) 



shown the possibility of representing each electronic state 
by a pair of classical action-angle variables. However, 
instead of these variables, cartesian electronic variables 
have been extensively used because of the simpler alge- 
braic form of the Hamiltonian. This approach and some 
of its extensions have been employed by several groups 
(for a complete list of references on the development and 
applications of this method, see [HI, [24| and references 
therein). Furthermore, by using the action-angle vari- 
ables, this formalism has also been recently implemented 
to study the dissipative dynamics (zero temperature) of 
chiral molecules. In this case, the asymmetry of the as- 
sumed double well potential model is due to the parity 
violating energy difference [251427} . A natural extension 
of this dynamics is to include the temperature effect along 
with a noise term. 

In this spirit, this is the first of a series of upcom- 
ing articles in which we describe, implement and apply 
a Langevin canonical approach to the dynamics of the 
TLS interacting with an environment at finite tempera- 
tures. Section II is devoted to introduce the canonical 
formalism for an isolated TLS, as well as its usual con- 
nection with the density matrix approach and the cor- 
responding thermodynamics. Finite temperature effects 
are included by means of noise-induced dynamics via a 
Caldcira-Lcggett-like Hamiltonian. Numerical simula- 
tions of this stochastic approach are presented and dis- 
cussed in Section III, where population differences (which 
are proportional to the optical activity) as well as co- 
herences are issued by assuming an Ohmic regime, for 
moderate-to-high temperatures. In particular, the com- 
petition between the tunneling and the asymmetry or 
bias displayed by the TLS is also analyzed. In all cases, 
the coherent and incoherent tunneling regimes are ana- 
lyzed and discussed in terms of NIBA path-integral ana- 
lytical expressions to provide the new parameters ruling 
this stochastic dynamics. Even more, the corresponding 
thermodynamic functions are obtained from this stochas- 
tic analysis at asymptotic times where the thermal equi- 
librium is reached. Finally, Section IV presents the main 
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conclusions and future perspectives of the formalism here 
employed. 



II. THEORY 

A. Canonical formalism for two level systems 

Let us consider an isolated TLS described by the 
Hamiltonian H — Sa x + ea z where a x ,z stand for the 
Pauli matrices. The isolated TLS is usually modelled 
in a phenomenological way by a two-well (asymmetric) 
potential within the Born-Oppenheimer approximation. 
From the knowledge of the eigenstates, |1), |2), the left 
and right states (or chiral states), \L) and \R), respec- 
tively, can be expressed by means of a rotation of angle 9 
given by tan26> = 6/e, where (L\H\R) = -6 (with S > 0) 
describes the tunneling rate and 2e = (L\H\L) — (R\H\R) 
(e can be positive or negative) accounts for the asymme- 
try due to the electroweak parity violation (for a chiral 
system) or any other bias term (for example, a magnetic 
field). 

Among other interesting representations of the isolated 
TLS [HI , an alternative and useful way of looking at it is 
based on the polar decomposition of the complex ampli- 
tudes entering the wave function (canonical description). 
The solutions of the time-dependent Schodinger equation 
{% = !) 



(1) 



can be written as |*(i)) = a L {t)\L) + a R (t)\R). If 
the complex amplitudes are written in polar form as 
a L,ii(t) = \aL,n(t)\e l ® L - R ( t \ and the population and 
phase differences between chiral states are defined as 
z(t) = \a R (t)\ 2 - \a L (t)\ 2 and $(t) = $ R (t) - $ L (t), 
respectively, it is an easy exercise to prove that the av- 
erage energy in the normalized | state is given by 
(V\H\y) = -25V1 - z 2 cos$ + 2ez = H , where H is a 
Hamiltonian function. As z and <3> can be seen as a pair of 
canonically conjugate variables, the Heisenbcrg equations 
of motion (which are formally identical to the Hamilton 
equations) are easily derived from z = —dH /d§ and 
$ = dHo/dz. Explicitly, the coupled equations describ- 
ing an isolated TLS in these canonical variables are 



= -26\/l - z 2 sin$ 



$ = 26- 



;Cos$ + 2e. 



(2) 



Thus, Eqs. ((2|) are totally equivalent to the usual time 
-dependent Schrodinger equation (p}. In fact, the exact 
quantum beating expression can be obtained by noting 
that Hq is a conserved magnitude (25j . For simplicity, 
the adimensional time t — >■ 2<5 i will be used in the rest of 
this work. This re-scaling implies that the Hamiltonian 
function Hq can be expressed again as 



Notice that the first term of the Hamiltonian function J3]) 
accounts for the tunneling process and the second one for 
the underlying asymmetry (due to a bias or the parity- 
violating energy difference) putting on evidence the two 
competing processes in this simple dynamics. 

Finally, the connection between the canonical and the 
density matrix formalism can be established as follows. 
Let p be the density matrix for the TLS, whose ma- 
trix elements are given by p RjR = \a R \ 2 , p L x = \a L \ 2 , 
Pl,r = a LO-* R and p R ,L = o, R a* L . On the other hand, 
from (*(i)|if|*(i)) = Ti:(pH) = H , one obtains that 
the time average values of the Pauli operators are given 
by 

(o-z)t = Pr,r ~ Pl.l = z 

(<7 X )t = PR.L + PL.R = ~yl - Z 2 COS$ 

(& y )t = ipR,L - ipL,R = yl - z 2 sin $, (4) 
which is consistent with (H) = S(a x ) + t{d z ) and 



(5) 



The time population difference can also be split into 
two components, a symmetric and antisymmetric one, 
under the inversion operation consisting of replacing e 
by -e. 



B. Stochastic dynamics 

The dynamics of the isolated TLS can be reduced to 
simply solve Eqs. ([2]) and then form appropriate combi- 
nations of z and $ to recover the populations and coher- 
ences. When dealing with interactions with the environ- 
ment consisting of a high number of degrees of freedom, 
more sophisticated theoretical approaches are needed. 
They can be widely classified into two frameworks ac- 
cording to the picture of quantum mechanics used: the 
density operator formalism (Schrodinger and/or inter- 
action picture) and the generalized Langevin equation 
(Heisenberg picture) [3, EH . Due to the fact the formal- 
ism used here to describe a TLS is very much attached to 
the definition of a Hamiltonian function ([3]), the second 
choice is much more convenient, apart from being much 
less employed in the theory of open quantum systems. 
Within this canonical formalism, a Caldeira-Leggett-like 
Hamiltonian, Q where a bilinear coupling between the 
TLS and the environment is assumed, is widely accepted 
in the literature. In particular, we have recently devel- 
oped this formalism to study the dissipativc dynamics of 
chiral systems [25M27| . 



H = - \/i - z 2 cos $ 



(3) 



As previously stated [25| , noting that 4> and z play the 
role of a generalized coordinate and momentum, respec- 
tively, one can introduce interactions with the environ- 
ment by means of a system-bath bilinear coupling via a 
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Caldcira-Lcggett-likc Hamiltonian expressed as 



regime, Eqs. ([7J) read now 



2 2 



Qf^CiXi + $ 2 ^]c 2 A 4 , 



(6) 



where the sums run over the coordinates of the bath os- 
cillators {pi, Xi} and Aj, Ci and Wj are suitable dimension- 
less constants representing generalized masses, couplings 
with the environment, and oscillator frequencies, respec- 
tively. Although the requirement of a bilinear coupling 
has been relaxed in a previous publication [27j . in this 
work we will retain it for simplicity. 

It should also be noticed that the usual spin-boson 
Hamiltonian has been implemented in the Meycr-Millcr- 
Stock-Thoss representation by coupling the bath posi- 
tion coordinate with the population difference of the TLS 
(see Eq. (2.11) of HH). On the contrary, in the approach 
here employed, it is the phase difference of the TLS the 
canonical variable which is coupled to the bath position 
coordinate, not the population difference. In any case, it 
is expected that thermal equilibrium is reached indepen- 
dent on the coupling scheme assumed. 

Within this scheme, the corresponding coupled 
Langevin-type dynamical equations are given by 



■ yl — z 2 sin $ 
j(t-t')$(t')$(t')dt' + £(*) 



i z e 

$ = ; COS$ + -. 

vT^ 2 5 



(7) 



where the time-dependent friction (damping kernel) is 
expressed as 



j(t) = ^A lC 2 coscJi(i-i') 



(8) 



and the fluctuation force or noise is given by 



£(t) =J2i c ^ (xi(0) cos u)it + pi (0) sin w^) 

- c 2 A 4 $(0)cosu; 4 t (9) 

which depends on the initial conditions of both the sys- 
tem and the bath. Taking the bath oscillators as classi- 
cal variables (classical noise) , moderate-to-high temper- 
ature regimes are expected to be properly described by 
this approach, as will be discussed through the rest of 
the manuscript. 

If a Markovian regime is assumed, the standard prop- 
erties of the fluctuation force (Gaussian white noise) are 
given by: (£(t))p = (zero average) and (£(0)£(f))/3 = 
mkBT^S(t) (delta-correlated) where /? = (fc^T) -1 , ks 
being Boltzmann's constant. The friction is then de- 
scribed by 7(t) = 2jS(t), where 7 is a constant and 
5(t) is Dirac's J-function (not to be confused with the 
(5-paramctcr describing the tunneling rate). Thus, in this 



= -v 7 ! - 2 2 sin$ - 7$(t) 



(10) 



The corresponding solutions provide stochastic trajecto- 
ries for the population and phase differences encoding 
all the information on the dynamics of the non-isolated 
TLS. These solutions are dependent on the four dimen- 
sional parameter space (e, <J, 7, T), apart from the initial 
conditions z(0) = z and $(0) = $o- In terms of the 
averages of Pauli operators, we have now the condition 



r v ? t 



(11) 



instead of condition (J5J. 

The use of classical noise imposes some restrictions on 
the range of temperatures where this approach remains 
valid. At high temperatures, 3> Tvy (or 7 _1 3> h/3) 
thermal effects are going to be predominant over quan- 
tum effects which become relevant, in general, at times 
of the order of or less than h/3, sometimes also called 
thermal time. However, in general, at very low tempera- 
tures, <C Tvy (or 7 _1 <C fr(3) the noise is colored and 
its correlation function is complex and our approach is 
no longer valid. The dissipative dynamics in the classical 
noise regime is obtained at zero noise or zero temperature 

3C I . This regime has already been considered elsewhere 

27 j. 

Averages of population and phase differences as a func- 
tion of time, (z(t)) and ($(f)) respectively, are issued 
from running a high number of stochastic trajectories. 
The role of initial conditions has been extensively dis- 
cussed in the literature (see, for example, 0,[2!|) and we 
are not going to deal with this important issue. In our 
dynamical study, the system will be prepared in one of 
the chiral states, left or right (z = 0.999 or —0.999 in 
order to avoid initial singularities), and the initial phase 
difference $0 will be uniformly distributed around the 
interval [— 2w, 2n]. An important test needed to validate 
the formulation and approach here employed is to re- 
cover the main equilibrium thermodynamics properties 
of the non-isolated TLS from the stochastic dynamics at 
asymptotic times. 



C. Thermodynamics of two-level systems 

A detailed analysis of the thermodynamics of non- 
interacting chiral molecules assuming a canonical dis- 
tribution was carried out elsewhere |3lj . In particu- 
lar, thermal averages of pseudoscalar operators were ex- 
tensively analyzed. The canonical thermal average of 
an observable X is defined as (X)p = Tr(ppX) where 
pp = Z~ x e~^ H °, H is given by © and Z is the par- 
tition function. For chiral states, Z = 2 cosh(/3A) with 
A = ^ S 2 + e 2 . Thus, the thermal average for the pop- 
ulation difference and coherences (in the L-R basis) are 
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easily calculated to be 

ipz)p = J tanh(/?A) 

(&*)p = ^tanh(/?A) (12) 

From the knowledge of the partition function, the re- 
maining equilibrium thcrmodynamical functions are eas- 
ily deduced such as the Helmholtz free energy, the en- 
tropy, the heat capacity, etc. From such analysis, a crit- 
ical temperature given by [3l[ 



is derived when z displays an inflection point and the heat 
capacity a maximum as a function of the temperature. At 
temperatures higher than T c , the effect of e is masked by 
thermal effects which tend to wash out the population 
difference z (racemization). At temperatures lower than 
T c , the value of the ratio e/S is critical. When this ra- 
tio is close to unity, z is determined by the competition 
between tunneling and asymmetry or bias. When it is 
much greater than one, the tunneling process plays a mi- 
nor role and z keeps more or less its initial value. Finally, 
when this ratio is much less than one, the racemization 
is always present. 

In a previous work [32j , it was showed that at very low 
temperatures, a chiral or two level system could display 
Bose-Einstein condensation as well as a discontinuity in 
the heat capacity (reduced temperatures fe^T/A < 1). 
In this work, we are going to assume that in this regime 
the Maxwell-Boltzman canonical distribution still ap- 
plies, instead of the Bose-Einstein statistics, in order to 
know if the thermodynamics can also be reached within 
this formalism with a classical noise. 

Finally, it is worth stressing that the thermodynamic 
functions are independent on the friction coefficient. 
Thus, our averages values of population differences issued 
from this stochastic dynamics is independent on the fric- 
tion coefficient as time goes to infinity (when the thermal 
equilibrium with the bath is reached). 

III. RESULTS 

A. Numerical details 

The general strategy consist of solving the pair of cou- 
pled equations (fT0| for the canonical variables under the 
action of a Gaussian white noise, which is implemented 
using an Ermak-like approach [33l . |34| . Notice that in 
the Langcvin like coupled equations to be solved, the 
noise term only appears in the equation of motion of 
the ^-variable. The phase variable has been taken to be 
uniformly distributed between —2n and 2ir. In order to 
avoid the singular behavior found at z — > ±1, the initial 
value for z will be chosen to be around ±0.999, very far 
from the equilibrium condition. Moreover, when running 



trajectories there are trajectories visiting " un-physical" 
regions where \z\ > 1. This drawback is mainly associ- 
ated with the intensity of the noise since, for large val- 
ues of it (which depends on both the temperature and 
the friction coefficient), the stochastic ^-trajectories be- 
come unbounded. To overcome this problem, we have im- 
plemented a reflecting condition such that, for instance, 
when the trajectory reaches z > 1. we change 1 its value 1 
to 2 — z. The time step used goes from 10~ 6 to 10 -4 (di- 
mensionless units) depending strongly on the regime to 
be explored in order to follow properly the corresponding 
dynamics: high or moderate temperature regimes and lo- 
calized or delocalized regime. As noted in [231, unstable 
trajectories can also be found for certain values of e, 5 
and 7 in the simple case of dissipative but non-noisy dy- 
namics. As this problem persists in case of dealing with 
stochastic trajectories, not every triple (e, S, 7) gives place 
to a stable trajectory. In these cases, the time evolution 
of individual trajectories is not possible and a previous 
stability analysis is mandatory. However, in the stable 
case, a satisfactory description of population differences 
and coherences has been achieved by running up to 10 4 
trajectories in all cases considered. 



B. Stochastic population differences and 
coherences 

In this subsection, thermal effects both in the popula- 
tion difference and coherences will be studied, extending 
previous results in which only dissipativ e dy namics (that 
is, at zero temperature) was considered [27(. Specifically, 
we will focus on the role of the critical temperature such 
that, roughly speaking, separates coherent from incoher- 
ent tunneling. Regarding the internal dynamics of the 
TLS, we start analyzing the delocalized regime where 
6 > e. In particular, we have taken 6 = 1 and e = 0.5. 
Given these values, the critical temperature is T c ~ 1 
in units of A. Thus, to study its effects in the dynam- 
ics, temperatures ranging from 200 to 0.4 (in units of A) 
have been considered. As the friction has been taken to 
be a constant value, 7 = 0.1, both moderate and high 
temperature regimes (fcsT ~ h'j and ksT 3> h'y, respec- 
tively) are covered. The propagation time step has been 
taken to be between 10 -6 (high temperatures) and 10 -4 
(moderate temperatures) . 

Trajectory averages (black dashed curves) of the pop- 
ulation differences, (z)(t), are plotted in Fig. (fTJ) for the 
range of temperatures previously given above. Left pan- 
els (a, c and e) show calculations for T = 200, 20 and 2, 
respectively Right panels (b, d and /) show the calcu- 
lations for T = 40,4 and 0.4. The system reaches the 
thermal equilibrium value, given by the first expression 
of Eq. (TT2"j) , at asymptotic times, which is plotted by a 
dotted line in all the panels. Thus, as expected, more 
time is needed to reach the thermodynamic equilibrium 
at lower temperatures since the coherent tunneling (the 
oscillation regime between the two states) is dominating 
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the dynamics. The incoherent tunneling which prevails 
at the two highest temperatures analyzed (panels a and 
b) leads to racemization very rapidly. When the tempcr- 
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Figure 1. Time dependence of (z) for different temperature 
values (in units of A): (a) 200, (b) 40, (c) 20, (d) 4, (e) 2 and 
(f) 0.4. Black dashed curves are issued from solving Eq. (|10p 
and red dashed curves from Eqs. l|14j> (incoherent regime) 
and (| 15p (coherent regime). 

ature approaches the critical temperature, T c , (panels d 
and e) the competition between the thermal and the in- 
ternal energy scale of the TLS (roughly given by fc^T and 
HA, respectively) gives place to the appearance of the os- 
cillating profile, fingerprint of coherent tunneling. When 
the temperature is again lowered (panel f), the internal 
structure dominates the dynamics, and more pronounced 
oscillations become more persistent in time. Thus, for a 
fixed value of the friction coefficient, the critical temper- 
ature roughly defines the coherent-incoherent transition. 

In order to extract more physical information from our 
approach, these stochastic calculations are fitted to NIBA 
analytical expressions describing both coherent and inco- 
herent tunneling. In the weak Ohmic damping limit and 
moderate-to-high temperature regime the path-integral 
results show a dependence on time according to Q 



"72 1 



(14) 



(z)(t) = (z)p + [l-(z) 

for the incoherent regime and 

(z)(t)=a 1 e-^ t + (z) p 

+ [(1 - ai - {z)p) cos Sit + a 2 sin Sit] e~ 7 *, (15) 

for the coherent regime. In these two equations, 
72, ai, <Z2, 7i, 7, CI and (z)p arc considered as free param- 
eters. Here {z)p is equivalent to the first expression of 



Eq. (jT2|) . Notice the fairly good fittings of the average of 
the stochastic ^-trajectories to Eqs. (fT4")) and (fl~5"j) (with 
dashed red line) for all panels of Fig. ((!} . SI can be inter- 
preted as the oscillation frequency of an effective damped 
harmonic oscillator and 71, 72 and 7 the different relax- 
ation rates in this very involved dynamics. In particular, 
72 gives the effective decay rate for the incoherent tun- 
neling and 71 and 7 two effective decays for the coherent 
tunneling with different weights. In other words, for this 
last regime, SI and 7 give us globally the two different 
time scales observed in this stochastic dynamics since 71 
gives the incoherent contribution in this coherent regime. 
In the NIBA framework, all of these parameters are re- 
lated but quite cumbersome expressions. 
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Figure 2. Time dependence of the coherence {S , x )p for differ- 
ent temperature values (in units of A): (a) 200, (b) 40, (c) 
20, (d) 4, (e) 2 and (f) 0.4. Black dashed curves are issued 
from solving Eq. (|10[) and red dashed curves from Eq. (I16[) . 



Let us consider now coherences. However, instead of 
analyzing the average values of the phase difference given 
by (3>)(i), we have calculated the average value of {a x ){t) 
issued from our stochastic trajectory analysis, the sec- 
ond expression of Eq. (|4]). These values are plotted in 
Fig. ([2]) for the same temperature values as reported in 
Fig. (|T|) (the same temperature value is assigned to each 
panel). Notice that the hat of the Pauli operator in the 
x-direction has been removed. Thus, we are considering 
{&x(t)) = (— \/l — z 2 cos $) and the stochastic behavior 
on time is given by the black dashed curves, as before. 
Similar comments with respect to the incoherent and co- 
herent regimes are also applicable here. Even more, in 
the weak Ohmic damping limit and modcrate-to-high 
temperature regime our stochastic values are again fitted 
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T 


e/S 


71 


n 


7 




4 


0.5 


0.0164 


-1.120 


0.015 




4 


0.5 


0.0162 


-1.120 


0.015 




4 


1.0 


0.0210 


-1.410 


0.013 


<*> 


4 


1.0 


0.0200 


-1.410 


0.012 




40 


1.5 


0.1220 


-1.800 


0.067 




40 


1.5 


0.1220 


-1.798 


0.072 




40 


1.0 


0.1060 


-1.407 


0.069 




40 


1.0 


0.1160 


-1.408 


0.076 



£= 1.5, 5= 1.0,7=0.01 



Table I. Numerical fitting of the average of stochastic trajec- 
tories corresponding to population differences and coherences 
(Eqs. (TT5} and (TTB)) ). For simplicity, only T = 4 and 40 and 
certain values of the ratio e/S are shown. All the quantities 
are dimensionless. 



to the expression provided by NIBA according to Q 

(cr x )(t) = b 1 e-^ t + (a x ) fj (16) 
+ [-(&i + (<t x )p) cos fit + b 2 sin tit] e -7 ', 

where b\, 62 and {&x)p & re again taken as free parameters 
and the last parameter corresponds to the second expres- 
sion of Eq. (H2J). The same interpretations of O, 71 and 
7 can be given as before. 

The consistency of the fittings in Figs. ([T]) and © is 
supported by the obtention of identical numerical values 
for those parameters which are common to both expres- 
sions in Eqs. (|15[) and ([TB")) . as shown in tabic Q] for only 
two temperatures, T = 4 and T = 40. The different time 
scales clearly manifest in such a table. 



C. Thermodynamics from dynamics 

Finally, it should be noticed that the formalism here 
presented describes accurately the time dependence of 
both the population differences and coherences since 
it is an alternative way to solve the time-dependent 
Schrodingcr equation. It is also worth noting the ca- 
pabilities of this approach to describe properly equilib- 
rium thermodynamics in a large interval of temperatures, 
ranging from moderate to high ones. In the very low- 
temperature regime, where fcgT <C hj, quantum noise 
effects are expected to occur. In this regime, the non- 
commutativity of the canonical variables describing the 
bath, [x u pi] ^ 0, would lead to [£(*),£(*')] + 0. Thus, 
within this range of temperatures, the approximation 
here employed should break up. However, in our study, 
we have assumed that the thermodynamical behavior at 
low temperatures is also described by a canonical distri- 
bution and correct values have also been obtained under 
this assumption. 

Having in mind this limitation, we note that clas- 
sical noise effects properly describe the region where 
both thermal and internal effects driving the dynamics 
of the TLS take place (in particular, for temperatures 




«T> =-^Tanh(SA) 

z eq A 




Figure 3. Thermal dependence of the population differences 
and coherences in a large range of temperatures. Stochastic 
calculations are shown with black points. The standard re- 
sults (Eq. (|12[) ~) are given with solid lines. The parameters 
of the corresponding dynamics are: e = 1.5, 5 — 1.0 and 
7 = 0.01. 



close to T c ). In addition, the effects of the localiza- 
tion/delocalization (given by the ratio between e and 5) 
are correctly taken into account, as can be seen in Figs. 
([3]) and In these two figures, the thermodynamical 
functions given by Eqs. (|12|) arc plotted (solid lines) for 
the two main regimes studied: e > S (Fig. (O), and 
e < S (Fig. ((H)). Similar results are reported for the 
regime e ~ S. The black points in both sets of figures are 
the time asymptotic values issued from the stochastic dy- 
namics. As can be seen, the agreement is fairly good at 
all temperatures studied. Thus, we can conclude that 
both the localized and dclocalizcd regimes are correctly 
described by the approach here presented. Clearly, the 
thermodynamical behavior is independent on the friction 
value chosen. 



IV. CONCLUSION 

In this work, a Langevin canonical approach has been 
developed to study the thermodynamics of two-level sys- 
tems. These canonical variables have been related with 
both the population differences and coherences of the sys- 
tem, showing that complete information about the dy- 
namics can be encoded in them. Thermal effects due 
to a noisy environment, modelled as a collection of clas- 
sical harmonic oscillators coupled linearly with the sys- 
tem via a Caldeira-Leggett-like Hamiltonian have been 
taken into account. Extensive tests in the propagation 
of both the population differences and coherences have 
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Figure 4. Thermal dependence of the population differences 
and coherences in a large range of temperatures. Stochastic 
calculations are shown with black points. The standard re- 
sults (Eq. (|12[) ) are given with solid lines. The parameters of 
the corresponding dynamics are: e = 0.1, 5 = 1.0 and 7 = 0.1. 



been carried out, showing that the implementation of 
the noise term in the dynamics is properly accounted 
for. This approach has allowed us to follow the evolution 
of the system towards the thermodynamic equilibrium 
with very good accuracy. Both coherent and incoherent 
tunneling regimes have been properly described by this 
formalism for an Ohmic environment in the moderate- 
to-high temperature regime. Moreover, these stochastic 
results have been fitted to analytical NIBA expressions to 
extract more physical information of the parameters rul- 
ing this dynamics. The regime of very low temperatures 
(much smaller than the critical temperature) which re- 
quires a generalized Langevin equation is not considered. 
An extension of this work including correlation functions 
and more thermodynamics functions such as heat capac- 
ity and entropy is currently in progress. 
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